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Abstract 

Within the framework of effective Lagrangians we calculate the free en- 
ergy density for an O(A^) antiferromagnet in 2+1 dimensions up to three- 
loop order in the perturbative expansion and derive the low-temperature 
series for various thermodynamic quantities. In particular, we show that 
the magnon-magnon interaction in the 0(3) antiferromagnet in d=2+l 
- the 0(3)-invariant quantum Heisenberg antiferromagnet on a square 
or a honeycomb lattice - is very weak and repulsive and manifests itself 
through a term proportional to five powers of the temperature in the free 
energy density. Remarkably, the corresponding coefficient is fully deter- 
mined by the leading-order effective Lagrangian jC^jj- and does not involve 
any higher order effective constants from C^j^j: related to the anisotropics 
of the lattice - the symmetries are thus very restrictive in d=2+l. We 
also compare our results that apply to O(A^) antiferromagnets in 2+1 di- 
mensions with the those for O(A^) antiferromagnets in 3+1 dimensions. 
The present work demonstrates the efficiency of the fully systematic effec- 
tive Lagrangian method in the condensed matter domain, which clearly 
proves to be superior to spin-wave theory. We would like to emphasize 
that the structure of the low-temperature series derived in the present 
work is model-independent and universal as it only relies on symmetry 
considerations. 
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1 Introduction 



With the present article we would like to further promote the effective Lagrangian 
method in the condensed matter domain - in particular, we would like to demon- 
strate its efficiency in describing the thermodynamic properties of systems exhibiting 
collective magnetic behavior. While the low-energy properties of antiferromagnets in 
dimension c?=3-|-l have previously been investigated within the effective Lagrangian 
framework [H, [il, 0, 0] > the present study we focus our attention on antiferromagnets 
in dimension d=2+l. A thorough analysis of these condensed matter systems, using 
effective field theory methods, was performed in Refs. However, in 

the present paper, we go one step further in the perturbative expansion, taking into 
account contributions to the free energy density up to three-loop order. Our calcula- 
tion applies to any system with a spontaneously broken internal symmetry 0{N) — > 
O(A^-l), provided that the corresponding leading-order effective Lagrangian can be 
brought to a Lorentz- invariant form. This system will be referred to as O(A^) antifer- 
romagnet. Throughout the paper, when talking about dimension, we always refer to 
the space-time dimension d = dg + 1, where dg is the spatial dimension. 

The attentive reader may particularly wonder why the quantum Heisenberg anti- 
ferromagnet, 

H = —j''^^ Sm ■ Sn , J = const. , (1.1) 

n.n. 

falls into the class of 0(3) antiferromagnets, i.e., represents a system described by a 
Lorentz-invariant leading-order effective Lagrangian. After all, the lattice structure of 
a solid singles out preferred directions, such that the effective Lagrangian in general 
is not even invariant under space rotations. In the case of a cubic lattice, however, 
the anisotropy only shows up at higher orders of the derivative expansion [8] - the 
discrete symmetries of the three-dimensional cubic crystal thus imply space rotation 
symmetry at leading order in the effective expansion. The same is true for an anti- 
ferromagnet defined on a square or a honeycomb lattice, which represents the system 
considered in this paper. Hence, the leading-order effective Lagrangian describing 
the quantum Heisenberg antiferromagnet on a cubic (d=3+l) or square/honeycomb 
(d=2-|-l) lattice is invariant under space rotations and can be brought to a (pseudo-) 
Lorentz-invariant form [2|: Antiferromagnetic spin- wave excitations exhibit relativis- 
tic kinematics, with the velocity of light replaced by the spin-wave velocity. As we 
will explain in detail later on, the spatial anisotropics which indeed start manifest- 
ing themselves at next-to-leading order in the effective Lagrangian, will not affect 
at all the main result of the present paper. Hence, a Lorentz-invariant framework, 
even at next-to-leading order of the derivative expansion, is perfectly justified in our 
calculation. 

Goldstone's theorem, which represents the basis of the systematic effective La- 
grangian method, states that, if the symmetry G = O(A^) of the Lagrangian is spon- 
taneously broken to H = O(iV-l), we must have A^-1 Goldstone bosons in the broken 
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phase {N > 2). For A^=3, these low-energy degrees of freedom are identified with 
the two spin-wave excitations - or the two independent magnon particles - in the 
spectrum of the 0(3) antiferromagnet. 

If the perturbations, which explicitly break the internal rotation symmetry O(A^) 
of the Lagrangian, are small, the corresponding Goldstone excitations remain light and 
dominate the low-energy behavior of the system. The effective Lagrangian method 
thus also applies to antiferromagnets in weak external fields. Goldstone's theorem 
guarantees that the Goldstone particles interact only weakly at low energies such that 
a systematic perturbative expansion in the momenta and the external fields can be 
performed. In the present work we perturbatively evaluate the partition function in 
a power series of the temperature, in order to obtain low-temperature theorems for 
various quantities of physical interest. 

We would like to emphasize that, from the effective Lagrangian perspective, the 
analysis of the low-energy properties of the system is approached from a unified and 
model-independent point of view, based on the spontaneously broken symmetry of the 
system. The method applies to any system where the Goldstone bosons are the only 
excitations without energy gap. The essential point is that the properties of these low- 
energy degrees of freedom and their mutual interaction are strongly constrained by the 
symmetries inherent in the underlying theory, such as the Heisenberg Hamiltonian - 
the specific nature of the underlying theory or model, however, is not important. For 
general pedagogic introductions to the effective Lagrangian technique see Refs. joj. 
Brief outlines of the method may be found in Refs. [IGl. For specific applications to 
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condensed matter systems the reader may consult Refs. [3|, 0, [8l,lll|, [l^, [131, [Ij, [151, 116 



Although our analysis is general, referring to any system which exhibits a spon- 
taneously broken symmetry 0{N) — >■ O(iV-l) and a Lorentz-invariant leading-order 
effective Lagrangian, our interest will primarily be devoted to the special case A^=3, 
which describes the 0(3) Heisenberg antiferromagnet in d=2+l defined on a square 
or a honeycomb lattice. Here the internal 0(3) spin symmetry of the isotropic Heisen- 
berg model is spontaneously broken by the ground state which displays a non-zero 
staggered magnetization. This system has been widely studied in condensed matter 
physics and it will be instructive to compare our results with the findings derived 
within the microscopic Heisenberg model. As we will see, the effective field theory 
approach is by far more powerful than spin-wave theory. 

Apart from our specific application concerning the partition function of the 0(3) 
antiferromagnet in d=2+l up to three-loop accuracy, we would like to point out that 
the effective field theory approach to this condensed matter system has also proven 
to be very efficient in other applications: In particular, in a recent publication on 
the constraint effective potential of the staggered magnetization of the 0(3) anti- 



ferromagnet [I5j, the quantitative correctness of the magnon effective field theory 
has been demonstrated in great detail at permille level accuracy by comparison with 
Monte Carlo simulations of the quantum Heisenberg model using the very efficient 
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loop-cluster algorithm. At this accuracy, three-loop effects clearly start manifesting 
themselves as there are small discrepancies between the Monte Carlo data and the 
two-loop predictions of the effective field theory. Indeed, it would be interesting to 
extend the finite-volume effective field theory formulas for the constraint effective 
potential to three loops and to confirm the correctness of the effective field theory ap- 
proach on an even higher level of accuracy. While this finite-volume calculation may 
be performed in a future study, in the present work we focus on finite temperature. 

The paper is organized as follows. For the sake of selfconsistency, in section [21 
we give a brief outline of the effective Lagrangian method at finite temperature. In 
section [3] we present the evaluation of the partition function up to three-loop order 
in the perturbative expansion. The issue of renormalization is then discussed in 
section m Section [5] contains our main results, i.e., the low-temperature series for the 
free energy density and other thermodynamic quantities up to three-loop order. In 
section [H] we justify why it is legitimate to use a (pseudo-)Lorentz-invariant framework 
in our calculation. We then compare in section [7| our results which apply to 0(3) 
antiferromagnets in three dimensions with those for 0(3) antiferromagnets in four 
dimensions. Finally, section [8] contains our conclusions, while some technical details 
concerning the renormalization and the evaluation of a specific three-loop graph are 
relegated to three appendices. 



2 The effective Lagrangian method at finite tem- 
perature 

In a Lorentz-invariant framework the construction of effective Lagrangians is straight- 
forward ITj: One writes down the most general expression consistent with Lorentz 
symmetry and the internal, spontaneously broken symmetry G of the underlying 
model in terms of Goldstone fields U"'{x),a = 1,..., dim(G)-dim(H) - the effec- 
tive Lagrangian then consists of a string of terms involving an increasing number of 
derivatives or, equivalently, amounts to an expansion in powers of the momentum. 
Furthermore, the effective Lagrangian method allows to systematically take into ac- 
count interactions which explicitly break the symmetry G of the underlying model, 
provided that they can be treated as perturbations. 

In the particular case we are considering, the symmetry G = O(A^) is explicitly 
broken by an external field. It is convenient to collect the {N-1) Goldstone fields U"" 
in a A^-dimensional vector W = {U^ ^U"") of unit length, 

U\x)U\x) = 1, (2.1) 

and to take the constant external field along the zeroth axis, if* = [H, 0, . . . , 0). The 
Euclidean form of the effective Lagrangian up to and including order p*^ then reads 
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Ill: 

-h^^iHWf - h^^WH\ (2.2) 

In the momentum power counting scheme, the field U{x) counts as a quantity of 
order 1. Derivatives correspond to one power of the momentum, oc p, whereas 
the external field H counts as a term of order p^. Hence, at leading order (oc p^) we 
have two coupling constants, F and E^, while at next-to-leading order (oc p^) already 
five constants, ci, 62, fci, and k^, show up. Note that these couplings are not fixed 
by symmetry - they parametrize the physics of the underlying theory and have to 
be determined either experimentally or in a numerical simulation. Using magnetic 
terminology, the square of the effective coupling constant F is the spin stiffness, 
while for the 0(3) antiferromagnet the quantities and W represent the staggered 
magnetization and the staggered external field, respectively. 

The effective Lagrangian method provides us with a simultaneous expansion of 
physical quantities in powers of the momenta and of the external field. The essential 
point is that, to a given order in the low-energy expansion, only a finite number of 
coupling constants and only a finite number of graphs contribute. The leading terms 
stem from tree graphs, whereas loop graphs only manifest themselves at higher orders 



in the derivative expansion 18 



A crucial difference with respect to the effective analysis in four space-time di- 
mensions concerns the suppression of loops in Feynman graphs: While loops are sup- 
pressed by two momentum powers in four space-time dimensions, in three space-time 
dimensions loop corrections are suppressed by only one power of momentum [l^ . As 
a consequence, the number of Feynman graphs that contribute to the perturbative 
expansion of the partition function up to a given order p", will depend on the space- 
time dimension. As we will see in the next section, there are fewer graphs in three 
dimensions that contribute up to three-loop order. 

The effective Lagrangian technique can readily be extended to finite temperature 



20|. In the partition function, contributions of massive particles are suppressed ex- 
ponentially, such that the Goldstone bosons dominate the properties of the system at 
low temperatures. In the power counting rules, the role of the external momenta is 
taken over by the temperature, which is treated as a small quantity of order p. The 
interaction among the Goldstone degrees of freedom in three dimensions generates 
corrections of order p/F cc T/F, while in four dimensions, the corrections are of order 
pV^' oc TyF\ 

In the effective Lagrangian framework at finite temperature, the partition function 
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is represented ClS db Euclidean functional integral 
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Tr [exp{-n/T)] = J [dU] exp J d^xC^ff) ■ (2.3) 

The integration is performed over all field configurations which are periodic in the Eu- 
clidean time direction, U{x^Xi + (3) = U{x, X4), with /5 = 1/T. The low-temperature 
expansion of the partition function is obtained by considering the fluctuations of the 
field U around the ground state V = (1, 0, . . . , 0), i.e. by expanding f/° in powers of 
U"", = a/1 — U"-U°-. The leading contribution (order p^) contains a term quadratic 
in U"" which describes free (pseudo-)Goldstone bosons of mass 

= T.sH/F'^ . (2.4) 

The remainder of the effective Lagrangian is treated as a perturbation. Evaluating 
the Gaussian integrals in the standard manner, one arrives at a set of Feynman rules 
which differ from the conventional rules of the effective Lagrangian method only in 
one respect: the periodicity condition imposed on the Goldstone field modifies the 
propagator. At finite temperature, the propagator is given by 

00 

G{x) = A{x,Xi + n(3), (2.5) 

n = —00 

where A{x) is the Euclidean propagator at zero temperature. We restrict ourselves 
to the infinite volume limit and evaluate the free energy density z, defined by 

z = -T lim L"Mn[Trexp(-7^/r)] . (2.6) 

L— >oo 

To evaluate the graphs of the effective theory, it is convenient to use dimensional 
regularization, since the symmetries of the theory are preserved within this scheme. 
The zero-temperature propagator then reads 

//•oo 
d'^pe'P"(M2 + p2)-i = y dp(47^p)-^/2e-^^^'-"'/^^ (2.7) 



3 Feynman graphs 

Our aim is to evaluate the partition function of an O(A^) antiferromagnet in dimension 
d=2+l up to three- loop order - the relevant Feynman graphs are shown in FiglH 

At leading order (order p^), we have a tree graph involving Clj-j-. The next order 
is p^, where we have a one- loop graph. Remember that in three dimensions every loop 
leads to a suppression of one momentum power only. At order p*^ the next-to-leading 
order Lagrangian C^rr contributes to a tree graph, while the leading Lagrangian C^rr 
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Figure 1: Feynman graphs related to the low-temperature expansion of the partition 
function for an O(A^) antiferromagnet up to three- loop order in dimension d=2+l. 
The numbers attached to the vertices refer to the piece of the effective Lagrangian 
they come from. Vertices associated with the leading term C^^j- are denoted by a dot. 
Note that loops are suppressed by one momentum power in d=2+l. 

manifests itself in the form of a two-loop graph. The situation is more involved at 
order p^, where we have three three- loop graphs with insertions from C^jj-, as well 
as a one-loop graph involving C^jj-. Note that higher-order pieces of the effective 
Lagrangian, starting with C^jf, are not relevant for the evaluation of the partition 
function in d=2+l at the three-loop level. 

In order to compare our three-loop calculation in 2-1-1 dimensions with the one 
referring to the evaluation of the partition function for an O(A^) antiferromagnet in 
3-1-1 dimensions, we have displayed the relevant graphs for the latter case in FigJ21 
Note that loops are now suppressed by two momentum powers, which leads to fourteen 
diagrams up to the three-loop level, whereas in three dimensions we only have eight 
diagrams. As we will see later on, the different loop counting and thus different 
organization of Feynman graphs reflects itself also in the renormalization of these 
graphs which turns out to be less involved in three dimensions. 

Let us now consider the case d=2+l and address the contributions from the rele- 
vant Feynman graphs of Fig{T] individually. At order p"^, the tree graph involving the 
leading-order effective Lagrangian leads to a temperature-independent contribution. 



At order p^, the one-loop graph involving C^jj-, which represents the free Bose gas 
term, is given by 



Z2 = -F^M^ . 



(3.1) 



Z3 - -o 



i(Ar-l)(47r)-^/'r(-|)M3 - 1{N - 1) go{M,T) . 



(3.2) 



7 



2 4A 4B 6A 6B 6C 




8A 8B 8C 8D 




8E 8F 8G 8H 



Figure 2: Feynman graphs related to the low-temperature expansion of the partition 
function for an 0{N) antiferromagnet up to three- loop order in dimension d=3+l. 
The numbers attached to the vertices refer to the piece of the effective Lagrangian 
they come from. Vertices associated with the leading term Cljrj- are denoted by a dot. 
Note that loops are suppressed by two momentum powers in d=3+l. Note also that 
we have used capital letters A, B, H in the definition of the diagrams, in order to 
distinguish them from the diagrams in d=2+l displayed in FiglH where we have used 
lower-case letters. 

The function go{M,T) is part of a set of kinematical functions go{M,T), gi{M,T) 
and g2{M,T) which are associated with the c?-dimensional noninteracting Bose gas 
and are defined by 

griM, T) = 2 / eM-pM') eM-n'/^pT') . (3.3) 

Jo (47rp) 

There are two graphs at order p^. The tree graph 4a involves two coupling con- 
stants from the next-to-leading order Lagrangian, 

z^a = -{k2 + h)M\ (3.4) 

while the two-loop graph 4b leads to a temperature-dependent contribution 

z,, = UN-l){N-3)^{G^)\ (3.5) 
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The expression Gi denotes the value of the thermal propagator at the origin 

Gi = G(x)|,=o- (3.6) 



The situation is more complicated at order p^, where we have four graphs: The 
three-loop graphs 5a, 5b, and 5c, involving the leading-order effective Lagrangian Cl^f, 
as well as a one-loop graph with a vertex from the next-to-leading order Lagrangian 

^eff 

Graph 5a factorizes into a term which only involves the thermal propagator at the 
origin: 
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z,a = UN +1){N- 1){N - 5)-^ . (3.7) 

Likewise, graph 5b exclusively contains propagators or derivatives thereof evalu- 
ated at the origin. 



z,, = -\{N -1){N- 3)^ {G,f - ^(7V -1){N- 2,f^ {G,f G, . (3.8) 



The quantity G2 corresponds to an integral over the torus T = 7?.^* x S^, with circle 
defined by -P/2<X4< P/2, and reads 



J d'^x^^G{x)y . (3.9) 



This integral can be expressed in terms of the derivative of the propagator at the 
origin with respect to the mass. 

Graph 5c leads to integrals over products of four propagators. Integrating by 
parts, they can be brought to the form 



= UN-l)iN-3)— J, - l(iV-l)(iV-2)— J, + iiV(iV-l)— {G,f , (3.11) 
where the functions Ji and J2 are given by 

Jl = j^d'^x [G{x)y , 

J2 = J^d''x[d^G{x)d^G{x)y . (3.12) 

Finally, the one-loop graph 5d with an insertion from ^^ff yields 

Z5d^{N-l){k2-h)^G,. (3.13) 
9 



Collecting all the pieces, we obtain the following expression for the free energy 
density of an O(A^) antiferromagnet in dimension d=2+l up to and including three 
loops: 

z = -F^M^ - |(iV-l)(47r)-'/'r(-|)M=^-i(iV-l)(7o(M,T) - (A;2 + fc3)M^ 
+ UN - l)(iV - 3)— {G^f + MN -1){N- 3)(3iV - 7)— {G^f 



- UN-l){N-?,f—{G^fG, + ^(iV-l)(iV-3)— Ji 



- i(iV-l)(iV-2)-lj2 + (iV-l)(A:2-fci)^Gi + 0(/). (3.14) 

Note that the quantities above involve the bare mass M of the (pseudo-)Goldstone 
bosons given in eq. fl2.4l) . The thermodynamics of the antiferromagnet is contained 
in the functions go,Gi,G2, Ji, and J2 which depend in a non-trivial manner on the 
ratio M/T. In the following section, we will take care of the singularities contained 
in the above expression and derive the low-temperature expansion for the free energy 
density. 



4 Divergences at d=2+l and renormalization 

In order to analyze the divergences in the limit c? — > 3, we split the thermal propagator 
into two pieces, 

G{x) = A{x) + G{x) , (4.1) 
where A(x) represents the propagator at zero temperature. At the origin, we have 

Gi = 2M'X + g,{M,T), 

G2 = (2-d)A + (?2(M,T). (4.2) 

The temperature-dependent quantities gr{M, T), defined in eq. (l3.3p . are smooth func- 
tions in the limit d 3. The temperature- independent contributions involve the 
parameter A 

A = i(47r)-''/'r(l-f)M''-^ (4.3) 
Remarkably, A is finite in the limit d 3, 

A = ^ • (4.4) 

On the other hand, in the limit (i — >■ 4 the parameter A contains a pole due to the 
singular behavior of the Gamma function. Accordingly, logarithmic divergences in 
the ultraviolet show up in four space-time dimensions, such that the next-to-leading 
order effective constants will undergo a logarithmic renormalization in d=3+l. For 
the moment, however, we focus on the case d=2+l. 
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In order to remove the singularities in the remaining integrals Ji and J2, as we 
show in appendix A, it suffices to subtract counterterms of the form ci + C2 gi{M, T) 
and C3 + C4gi{M,T), respectively, 

Ji = Ji - ci - C2fi'i(M,T) , 

J2 = J2-C3-C4(7i(M,T), (4.5) 

where the constants Cj are singular functions of the dimension d. While the quantities 
Ci and C3 renormalize the vacuum energy, C2 and C4 renormalize the mass M (see 
below). 

We now insert the decompositions (14. 2 p and (14. 5 p into the free energy density 
(I3.14P and discuss the various pieces therein. All contributions which are independent 
of the temperature, 

zo = -F'M'-^{N-l)M'-{h + h)M' + ^{N-l){N-3) — 
- -1){N- 3){9N - 23)^ + ^(AT - 1){N - 3)^ ci 



1 



- liN -l){N-2)—cs-liN~ -h)— + 0{p') , (4.6) 



merely renormalize the vacuum energy. 

Next, we consider all terms in the free energy density (I3.14p which are linear in the 
kinematical functions gr{M,T). In appendix B we show that these contributions can 
be merged into a single term proportional to qqIMt^jT) by renormalizing the mass, 
M —>■ M^, according to 

= + (AT - 3) A 1; + {2(i. - + I5 + '-^} ^ + O(M') . (4.7) 

The quantity bi is related to the singularities contained in the coefficients C2 and C4. 
We thus see that the divergences in 61, originating from the three- loop graph 5c, are 
absorbed into the combination /c2 (stemming from the one-loop graph 5d) of next- 
to-leading order effective constants. After mass renormalization, the only surviving 
term linear in the kinematical function is the contribution from the free energy density 
of noninteracting magnons given by 

-l{N-l)go{M^,T), (4.8) 

which now depends on the renormalized mass M,r- 

Finally we have to take care of the terms quadratic and cubic in the functions 
gr{M,T), which is also done in appendix B. We are then left with the following 
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expression for the free energy density of an O(A^) antiferromagnet in three dimensions, 



z 



_^(iV - l)(iV - 3)(5iV - 11) {g,f + ^(.v - x^.v - ono.v -nyiKyi 



^(iV - l)(iV - 3)(5iV - 11) ((7i)' + - l)(iV - 3)(3iV - 7) ((7 

- UN-l){N-3f^{g,fg, + ^^+0{p'), (4.9) 

where we have defined the function Q{Mt^, T) by 

Q^i^iN- 1){N - 3) Mt Ji - - l)(iV - 2) J2 . (4.10) 

The expression (14. 9 p for the free energy density is free of divergences and only involves 
the physical mass M^. In particular, the kinematical functions are defined as gr = 

For dimensional reasons, the thermodynamic functions in eq. fl4.9p are of the form 
T^f^r), where r is the dimensionless ratio 



T 

m: 



(4.11) 



Explicitly, in d=2+l they are given by 



go = T'ho{T), g, = Th,{r), g, = ^h2{T), Q = q{T) , (4.12) 



such that the free energy density can be written as 



1 



- T2b(^-l)(^-3)(5iV-ll)Tij/^iW'2^' 



1 



- UN-l){N-3f^^h{Tfh,{T)r^ + j-^q{T)T' + 0{T'). (4.13) 

This expression for the free energy density of an O(A^) antiferromagnet in 2+1 dimen- 
sions represents the basic result of our paper. The ratio r = T/M^r can take any value, 
as long as the quantities T and M^^ themselves are small compared to the intrinsic 
scale A of the theory which, in the case of the 0(3) antiferromagnet, may be identified 
with the exchange integral J of the Heisenberg model fll.ip . 

Remarkably, for = 3 - the quantum Heisenberg antiferromagnet on a square 
lattice - most of the terms drop out and we are left with the following simple expression 
for the free energy density of the 0(3) antiferromagnet in d=2+l: 

z = zo- hoir) T' + ^ g(r) + ^(T^) (iV = 3) . (4.14) 

12 



While the term cubic in the temperature corresponds to the free Bose gas, the term 
proportional to five powers of the temperature represents the leading contribution due 
to the spin-wave interaction. Note that for this special case (A^=3), the function g(r) 
defined in eq.f l4.12p only involves the contribution proportional to J2. 



5 Low-temperature series for the 0(3) antiferro- 
magnet in d=2+l 

With the representation fl4.14p for the free energy density of the 0(3) antiferromagnet 
in 2+1 dimensions, we are now able to discuss various thermodynamic quantities for 
this system. We are particularly interested in the limit T 3> M^^ which we implement 
by holding T fixed and sending M^^ (or, equivalently, the external field H) to zero. 
Since we keep the fixed T small compared to the intrinsic scale A of the underlying 
theory, we do not leave the domain of validity of the low-temperature expansion. 

Because the system is homogeneous, the pressure is given by the temperature- 
dependent part of the free energy density, 

P = Zo-z = hoir) T3 - i- g(r) T' + 0{T') . (5.1) 



The non-trivial dependence of the quantity P on the ratio r = T /M^^ is contained 
in the functions ho^Mj,, T) and q^Mt,, T), which are defined in eqs.f l3.3p and fl4.12p . For 
the function /io(Mn-, T) an analytical expression can be provided in the limit T ^ 
(see appendix C): 



T) = - 

TT 



1M2 1M2 



1M3 1 Mt 

TT I TV 

6 T3" 96 T4 



+ O 



T 



(5.2) 



The function q{Mj^,T), on the other hand, we have evaluated numerically, using 
the representation for J2 given in appendix A - a plot of is provided in Figj3l Still, 
in the limit T 3> M^, the function may be parametrized by 

,«^...#.o(^)\ .^-i-, (.3, 

where the coefficients qi are real numbers. Here we only need qi which takes the value 

= -0.008779 . (5.4) 
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Figure 3: The function q{a) for A^=3, where a is the dimensionless parameter a 
M^/2ttT = l/27rr, introduced in Ref. 



Making use of the above representations for HqIMt^jT) and q{MT^,T), in the hmit 
T ^ the pressure takes the form 



P 



03)^3 



TT 



0.3826 



1 + 0.02294— + 0(r3) . (5.5) 
F 



The corresponding expressions for the energy density u, for the entropy density 
s, and for the heat capacity cy for the 0(3) antiferromagnet in 2+1 dimensions, are 
readily worked out from the thermodynamic relations 



5 u = Ts — F . 
dT ' 



cv 



du ds 
df ^ df 



(5.6) 
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with the result 



u 



cv 



203)^3 

TT 

0.7653 
303) t2 

TT 

1.1479 
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TT 

2.2958 



C(3) ^ ' 

1 + 0.04589 — + 0{T^) 



1 



303) ^ 



1 + 0.03824— + 0(r3) 



3C(3) F4 
1 + 0.07648 — + 0(r3 



(5.7) 



The respective first terms in the above series represent the free Bose gas contribu- 
tion which originates from a one-loop graph. The effective interaction among the 
Goldstone bosons only manifests itself through a term of order in the pressure, 
related to a three-loop graph. Interestingly, the coefficient qi is negative, such that the 
magnon-magnon-interaction in the 0(3) antiferromagnet in d=2+l is repulsive at low 
temperatures. It is remarkable that the coefficient of the interaction term in these 
series is fully determined by the symmetries inherent in the leading-order effective 
Lagrangian, and does not involve any next-to-leading order coupling constants from 
jO-^ff, reflecting the anisotropics of the square lattice or the Lorentz-noninvariant na- 
ture of the quantum Heisenberg antiferromagnet defined on a square or a honeycomb 
lattice - the symmetry is thus very restrictive in d=2+l. Note that there is no inter- 
action term of order T*^ in the pressure: the two-loop contribution Z4b is proportional 
to — 3 and thus vanishes for the 0(3) antiferromagnet, irrespective of the actual 
value of the ratio r = T/M^.. 

The fact that an interaction term proportional to four powers of the temperature 
does not show up in the temperature expansion for the pressure of the 0(3) antiferro- 
magnet in the limit T ^ Mj^was already pointed out in Ref. 0]: this was an effective 
Lagrangian calculation that operated on the two-loop level. We are not aware of any 
microscopic calculation that aimed at this accuracy. Moreover, our result that the 
leading contribution of the magnon-magnon interaction in the pressure is repulsive 
and of order requires a three-loop calculation on the effective level, performed in 
the present study - it is probably fair to say that this accuracy is beyond the reach of 
any realistic microscopic calculation based on spin-wave theory. 

In fact, as pointed out in Ref. [8|, there were inconsistencies between the results 
obtained by spin-wave theory, Schwinger boson mean field theory and Monte Carlo 
simulations, already with respect to the very leading term (oc T^) in the temperature 



expansion of the energy density 2J, |25|, l26|] . The error was later attributed to some 



numerical problems in solving the equations arising in Schwinger boson mean field 
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theory. The systematic effective Lagrangian method, which approaches the problem 
from a unified and model-independent perspective based on the symmetries of the 
underlying theory, thus clearly proves to be superior to these conventional condensed 
matter techniques. 

As demonstrated in a recent publication on the constraint effective potential of the 



staggered magnetization of the 0(3) antiferromagnet [15|], three-loop effects clearly 
start manifesting themselves, as there are small discrepancies between the very precise 
Monte Carlo data and the two-loop predictions of the effective field theory. Indeed, 
it would be interesting to extend the finite-volume effective field theory formulas for 
the constraint effective potential to three loops, in order to confirm the correctness of 
the effective field theory approach on an even higher level of accuracy and to extract 
the numerical values of some combinations of effective next-to-leading order coupling 
constants. Although non-trivial, this would certainly be feasible within the framework 
of the systematic magnon effective field theory. 

We now make an important comment on the range of validity of the above low- 
temperature series. After all, we are considering the limit T ^ M^-, which we have 
implemented by holding T fixed and sending M^^, or, equivalently, the external field 
H, to zero. However, following Mermin and Wagner |27(], there is no spontaneous 
symmetry breaking at any finite temperature in the 0(3)-invariant Heisenberg model. 
Accordingly, there are no massless magnons in the low-energy spectrum at any fi- 
nite temperature. Rather, the magnons pick up an exponentially small mass. The 
argument of the exponential is proportional to the inverse temperature, 

,n = (-) 2rf^ exp[ - — ] 1 - 5^ + O(^) . (5,8) 



as derived in Ref. [5|, 128[. Strictly speaking, it is therefore not legitimate to switch 
off the external field H completely, because the above calculation does not take into 
account the non-perturbative effect of m. However, the corrections due to the non- 
perturbatively generated mass gap are so tiny, that they cannot manifest themselves 
in the above power series. In order to verify this claim, we now estimate the order of 
magnitude of these corrections. 

The above low-temperature series are valid as long as the correlation length of the 
Goldstone bosons, ^ = l/M^- is much smaller than the non-perturbatively generated 
correlation length ^„p = 1/m, let's say. 



1 _ ^ _ 8 2nF^ T 
1000 ~ ^ ~ e T 



27cF 



2-, 



(5.9) 



where we have used eq.f l5.8p . For the Heisenberg antiferromagnet in d=2+l, the spin 



stiffness F"^ has been determined very precisely in Monte Carlo simulations 15|, |29 
In units of the exchange integral J it takes the value F"^ = 0.1808(4) J, such that the 
quantity 27rF^ is of the order of J. Now the exchange integral defines a scale in the 
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underlying theory and for the effective expansion to be consistent, the temperature 
has to be small with respect to this scale. Assuming that 



T 1 

— ^ = — , (5.10) 
27rF2 100 ' ^ ' 

relation (15.91) then yields the ratio 

^^10-38. (5.11) 

Remember that, in the above low-temperature series, we have implemented the limit 
T ^ Mt, by holding T fixed and sending to zero. We thus see that, in principle, 
we cannot completely switch off the mass M^^ - rather, we start running into trouble 
as soon as the ratio M^^/T is of the order of the above value. However, the error 
introduced is indeed very small: the leading one-loop contribution in the free energy 
density, according to eqs. (14.91) and (15. 2p is 



I I Ml 
^4 T2" 



1 1 ^■ 

rp2 



10-'% (5.12) 



which is extremely small also with respect to the three-loop contributions in eq. (l5.5p . 
We thus confirm that the corrections due to the non-perturbatively generated mass 
gap are so tiny that they cannot manifest themselves in the above low-temperature 
expansions for the thermodynamic quantities. In other words, the subtleties raised 
by the Mermin- Wagner theorem in d=2+l are not relevant for our calculation. 

While our effective calculation is restricted to the regime ^ ^ ,^„p, the regime 
i 3> ^np, is perfectly well accessible also with effective field theory methods. However, 
one has to resort to a different type of perturbative expansion. A similar situation 
occurs when one considers finite size effects: when the Goldstone boson mass is small 
compared to the inverse size of the box, a different effective expansion scheme, the 
so-called e-expansion, applies. Indeed, various problems within this framework have 



been investigated in detail [30|, l3l|, l32|, |33 . 



6 Justification of the Lorentz- invariant framework 

In this section we would like to explain why it is justified to use a Lorentz-invariant 
framework in our calculation. We have seen that anisotropics induced by a cubic or a 
square lattice do not affect the leading-order effective Lagrangian, such that Cl^^ can 
be written in a Lorentz-invariant form, where the spin-wave velocity takes over the 
role of the velocity of light: The accidental 0(3) space rotation symmetry at the Clff 
level implies (pseudo-) Lorentz invariance. On the other hand, the anisotropics related 
to the lattice structure do show up in the next-to-leading order effective Lagrangian 
£g r r. Now, in our three-loop calculation of the partition function in d=2+l the only 
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temperature-dependent diagram involving C^ff is the one-loop diagram 5d, which is 
quadratic in the magnon field. Indeed, here we have a new term due to the lattice 
anisotropies, also contributing to the diagram. The respective term, 

^ dsdsU'dsdsU\ (6.1) 

s=l,2 

is invariant under the 90 degrees spatial rotation symmetry of the square lattice, but 
not invariant under continuous 0(3) space rotations. Interestingly, this term is absent 
in the case of the honeycomb lattice, as it is not allowed by the 60 degrees rotation 
symmetry. 

However, both for the square and the honeycomb lattice there are additional terms 
showing up at next-to-leading order: If we consider an 0(3)-symmetric, i.e., space- 
rotation symmetric Lagrangian C^ff - and not a Lorentz-invariant Lagrangian C'^ff 
as we have done so far - there are further terms like 

AWAU' , WU'drU^drU'' , (6.2) 

that also have to be taken into account in C'^ff- The essential observation, however, is 
that all these Lorentz-noninvariant extra terms in eqs. fl6.ip and fl6.2p contributing to 
the one-loop graph 5d merely modify mass renormalization or give rise to higher-order 
corrections of the dispersion law, 

uik) =v\k\ + 0{P), (6.3) 

but cannot manifest themselves in the magnon-magnon interaction up to the order 
considered in the present work. Although they give rise to an additional term 
in the free energy density involving five powers of the temperature, this is a purely 
kinematical effect related to the one-loop graph 5d - the leading contribution due 
to the magnon-magnon interaction, also of order T^, will not be affected. Hence, 
our main result regarding the weakness and the repulsive character of the magnon- 
magnon interaction in the 0(3) antiferromagnet in d=2+l perfectly well applies to 
the quantum Heisenberg antiferromagnet defined on a square or a honeycomb lattice. 



7 Antiferromagnets in 2+1 and 3+1 dimensions: 
Structure of the low-temperature series 

In this section we want to compare the low-temperature series for antiferromagnets 
in 2+1 and 3+1 dimensions, pointing out differences as well as similarities. The 
effective Lagrangian method is ideally suited to understand the structure of these 
low-temperature series, as it adopts a unified perspective based on symmetry consid- 
erations. 
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As we have discussed in Sec|3l loops in four dimensions are suppressed by two 
momentum powers, whereas in three dimensions they are suppressed by one power 
of momentum only. Consequently, the organization of the loop expansion for the 
partition function depends on the space-time dimension and reflects itself also in 
the number of Feynman diagrams that have to be evaluated. Up to the three-loop 
level, we have the fourteen diagrams in four dimensions, displayed in Figj2]- in three 
dimensions there are only eight, displayed in FiglH In particular, in d=2+l, there are 
no two- loop graphs involving the next-to-leading order Lagrangian C^ff- Moreover, 
contributions from /^e// ^eff needed in d=2+l. One thus notices that the 

restrictions imposed by symmetry are extremely strong in d=2+l: Up to the three- 
loop level, no effective coupling constants from C^jj- or Clj^j enter the calculation and 
the couplings in - as we have seen - do not affect at all the spin-wave interaction 
part in the free energy density. 

Another immediate consequence of the dimension-dependent loop counting is the 
fact that interactions among antiferromagnetic magnons in d=2+l generate cor- 
rections of order p/F oc T/F, whereas in d=3+l these corrections are of order 
/ F"^ oc T'^/F'^. The low-temperature series for the various thermodynamic quanti- 
ties are thus expected to proceed in steps of one power of T in d=2+l and in steps 
of in d=?>+l - we will come back to this point below. 

We now briefly review the relevant results for an 0{N) antiferromagnet in 3-1-1 
dimensions - details of the calculation can be found in Ref. The formula for the 
pressure takes the form 



P = l{N-l)go + A7ra{giy + ng 



+ 0{p 



{d = 3 + l). (7.1) 



The temperature dependence is contained in the kinematical functions QriM^^^T) and 
in j(M^,T). In the limit H —>■ (or, equivalently, T ^ M^) we are interested in, 
analytical expressions for the functions go, gi and g can be provided (see Ref. [23| or 
appendix C), 



45 ^ 



1 - 



15 Ml 15 Ml 45(7 



+ 



47r2 T2 27r3 



+ 



ln47r)M^ 



167r4 T4 
(d = 3 + 1) 



J_rp2 



Ml, M^ 



675^ ^ 



47r2 T2 ^ \T 



(7.2) 



whereas the function j, containing the three-loop contribution from graph 8C, has to 



19 



be evaluated numerically, 



J = . lu- + + + 0{-f) . . . ^. (7,3) 

The coefficients ji and j2 in this expansion are real numbers. Note that the function 
j(r) diverges logarithmically in the limit H —* 0. It should be pointed out that the 
renormalized or physical mass M^^ in d=3+l is given by 



Ml 



( ^ M^ M^ 

M'+ [iN-3)X + 2ik2-h)}— + c— + OiM') (rf = 3 + l), (7.4) 



which is different from the analogous expression (14.71) in d=2+l, due to the different 
loop counting. Without going into details, we just mention that the quantity A, 
according to eg. (14.31) . is divergent in d=3+l and that the corresponding singularity 
occurring in the two- loop graph 6A (see Fig|2]) is absorbed into the combination ^2 — ^1 
of next-to-leading order coupling constants originating from graph 6B. The absorbtion 
of the divergences showing up in the various graphs of order p^, on the other hand, 
even involves a coupling constant of jC^jj stemming from graph 8G, which is contained 
in the quantity c of eg. (17.41) . 

Finally, the constants a and b in the pressure (17. ip involve the scales Ha and Hh, 
(iV-l)(iV-3)S,/7 {N-lf {i:,Hf ^^H 



32n F4 2567r3 ' 

which are related to coupling constants of C^jjr (for details see the appendix in Ref. [3]). 
The first term in a, linear in the external field H, originates from the two-loop graph 
6A. The logarithmic contributions in a and b involving the two scales Ha and Hb, 
originate from two-loop graphs with insertions from C'^jf- Note that in d=2+l these 
two-loop graphs are already beyond the next-to-next-to-leading order considered in 
the present paper. 

Eguipped with the above formulas, the low-temperature expansion of the pressure 
for an O(A^) antiferromagnet in d=3+l in the limit H amounts to 



90 



TV — 2 T 
72 T ^ ^ 



{d = 3 + l). (7.6) 



The first contribution represents the free Bose gas term which originates from a one- 
loop graph, whereas the effective interaction among the Goldstone bosons, remark- 
ably, only manifests itself through a term of order T®. This contribution contains 
a logarithm, characteristic of the effective Lagrangian method in four space-time di- 
mensions, which involves a scale, Tp, related to Hb (see the appendix in Ref. 0]). The 
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occurrence of a scale involving coupling constants from C^ff is a consequence of the 
space-time dimension d=3+l: In four dimensions the parameter A, defined in eq. fl4.3p . 
contains a pole, which can be absorbed into coupling constants of C^jj- by a suitable 
logarithmic renormalization. Note that the divergences in the function j eq. (17.31) and 
in the constant b eq. (17.51) cancel, such that the remaining expression involving the 
scale Tp is well defined in the limit H ^ 0. 

At low temperatures, the logarithm ln[Tp/T] in the pressure (17.61) is positive, such 
that the interaction among the Goldstone bosons in d=3+l, in the absence of an 
external field H, is repulsive, much like in d=2+l. The symmetries in d=3+l, how- 
ever, are somewhat less restrictive than in d=2+l, where the interaction term - the 
last term in eq. (l4.13l) involving the function g(r) - is unambiguously determined by 
the coupling constant F of the leading order effective Lagrangian: In d=3+l, next- 
to-leading order effective constants from C^jj- do show up in the scale Tp. Still, the 
symmetry is also rather restrictive in d=3+l, as it unambiguously fixes the coefficient 
in front of the logarithm in terms of the coupling constant F. Note that there is no 
term of order in the above series for the pressure. This is due to the fact that the 
respective two-loop contribution (graph 6A, Figl2l) in eq. (l7.1l) is proportional to the 
constant a that vanishes for a zero external field. 

Finally, the energy density u, the entropy density s, and the heat capacity cy in 
the limit H ^ are given by 
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cv 
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2 7r2(Ar-l)T=^ 
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1 + ^^f71n^-l) + 0(T6) 
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{d = 3 + l), 
(7.7) 



Note that the limit H can readily be taken in d=3+l, since the Mermin- Wagner 
theorem does not apply here: there are no exponentially small non-perturbative cor- 
rections in the above low-temperature series. On general power counting grounds one 
would expect the low-temperature series to proceed in steps of in 3-1-1 dimensions. 
However, as for the pressure before, there are no correction terms proportional to e.g. 
six powers of the temperature in the internal energy because the constant a vanishes 
in the limit H 0. 



For the specific case A^=3, the above series for the thermodynamic quantities take 
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the form 
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These series, vahd for the 0(3) antiferromagnet in (i=3+l, we now want to compare 
with the analogous series for d=2+l. 

As we have seen, in three dimensions the parameter A is finite and therefore no 
such scale, involving next-to-leading order coupling constants from C'^fj-, arises in the 
low-temperature series of the thermodynamic quantities. In the limit H ^ 0, the low- 
temperature expansions of the thermodynamic quantities for the 0(3) antiferromagnet 
in d=2+l take the form 
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(7.9) 



Here we would expect the low-temperature series to proceed in steps of T, since every 
loop in d=2+l leads to an additional suppression of one power of the temperature. 
However, as it was the case for d=3+l, there are no next-to-leading order corrections 
in the above series for d=2+l: a term proportional to T* in the pressure is absent. 
Again, the corresponding two-loop graph 4b is proportional to (see eq. (l3.5p ). such 
that it vanishes in the limit H ^ 0, much like the constant a in 3+1 dimensions before. 
In the absence of a staggered field, the magnon-magnon interaction both in d=2+l 
and d=3+l thus becomes very weak, as we are dealing with a next-to-next-to-leading 
order effect. 

For non-zero external field, the low-temperature representations of the thermody- 
namic quantities retain their form, except that the coefficients now become functions 
of Mtt/T. In the region T 3> one recovers the results of the theory for zero exter- 
nal field, whereas in the opposite limit, T <C M^r, the gas is dilute and the particles 
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move non-relativistically. The properties of the system are therefore very sensitive to 
the value of the ratio M^^/T. 

To illustrate this sensitivity, let us consider the pressure and discuss the general 
situation for N > 2. In the limit i7 — > 0, as we have seen, a two-loop contribution 
of order in d—S+1 - or in d—2+1 - does not occur. This is no longer the case 
for an approximate symmetry {H ^ 0): remarkably, the sign of the corresponding 
interaction term of order (oc HT^, graph 6A) in d=3+l - or the corresponding 
interaction term of order (oc HT^, graph 4b) in d=2+l - turns out to be negative. 
With respect to the limit i7 — > 0, the sign of this interaction term is thus different: 
in the absence of an external field, the first non-leading term (order p^ in d—2+1, 
order p^ in d—3+1) is positive and the interaction among the Goldstone bosons thus 
repulsive. We conclude that a weak external field damps this repulsion among the 
Goldstone bosons, such that the effective interaction becomes even weaker. 

Interestingly, the case A'"=3 is rather special: since the two-loop contribution both 
in d—2+1 and d=3+l is proportional to (A^-3), the above mentioned damping of the 
interaction does not occur. Still, the repulsive interaction between antiferromagnetic 
magnons in three or four dimensions is very weak as we are dealing with a next-to- 
next-to-leading order effect. 



8 Conclusions 

Condensed matter systems exhibiting a spontaneously broken continuous symmetry 
may very efficiently be analyzed with the fully systematic method of effective La- 
grangians. In the present study we have considered O(A^) antiferromagnets in d=2+l 
space-time dimensions which display a spontaneously broken internal rotation symme- 
try 0{N) — > O(A^-l) and whose leading-order effective Lagrangian can be brought to 
(pseudo-)Lorentz- invariant form. The low-temperature properties of this system are 
dominated by the corresponding Goldstone bosons, which for A'"=3 may be identified 
with the two antiferromagnetic magnons or spin-wave excitations. 

We have extended previous results for O(A^) antiferromagnets in d—2+1 to higher 
orders in the derivative expansion, evaluating the partition function up to and includ- 
ing three-loop diagrams. Although the renormalization and the subsequent numerical 
evaluation of one particular three-loop graph turns out to be non-trivial, the calcu- 
lation is perfectly feasible within the effective field theory framework. One of our 
main results is that the interaction among magnons in the 0(3) -invariant Heisenberg 
antiferromagnet, defined on a square or a honeycomb lattice, is very weak and re- 
pulsive at low temperatures, manifesting itself through a term proportional to five 
powers of the temperature in the pressure. Remarkably, the coefficient of this interac- 
tion term is fully determined by the leading-order effective Lagrangian Clj-j and does 
not involve any higher order effective constants from C,^tf - the symmetry is thus 
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very restrictive in d=2+l. As we have argued, additional effective constants in C^ff, 
talcing into account the Lorentz-noninvariant nature of the system, merely affect the 
renormalization of the magnon mass or yield higher-order corrections to the magnon 
dispersion law, but do not affect at all the leading contribution originating from the 
magnon-magnon interaction in the pressure. 

The free energy density for O(A^) antiferromagnets in d=3+l up to the three- loop 
level, on the other hand, does involve coupling constants from C^ff, which undergo 
logarithmic renormalization. Accordingly, the effective expansion of thermodynamic 
quantities now contains a scale, related to these coupling constants in £gjj. In four 
dimensions, within the (pseudo-)Lorentz-invariant framework, we thus need more phe- 
nomenological input, i.e. the numerical values of some next-to-leading order effective 
coupling constants, in order to fully specify the structure of the magnon-magnon in- 
teraction in the low-temperature series up to the three-loop level. Still, the symmetry 
is also very restrictive here, as it unambiguously fixes the coefficients in the expansion 
of the free energy density of 0{N) antiferromagnets in d=3+l up to order T^, where 
the logarithm involving the scale enters. 

The low-temperature theorems for the various thermodynamic quantities of 0{N) 
antiferromagnets in three and four space-time dimensions are exact up to and in- 
cluding three loops: independently of the specific underlying model, they are valid 
for any system with a spontaneously broken symmetry O(A^) — > O(A^-l), provided 
that the system can be described in a (pseudo-)Lorentz-invariant framework, with the 
velocity of light replaced by the spin- wave velocity. In particular, there are no approx- 
imations or idealizations involved in our main result regarding the weakness and the 
repulsive character of the magnon-magnon interaction in the 0(3) antiferromagnet 
in d=2+l: Although we use a (pseudo-)Lorentz- invariant framework, our calcula- 
tion is not just some kind of 'academic' exercise, as the lattice anisotropics, or the 
Lorentz-noninvariant nature of the system in general, cannot manifest themselves in 
the magnon-magnon interaction up to the three-loop order of the perturbative expan- 
sion, considered in the present paper - hence our calculation applies, as it stands, to 
the quantum Heisenberg antiferromagnet defined on a square or a honeycomb lattice. 

We would like to emphasize that the order of the calculation presented here, as we 
have argued in section [5], appears to be beyond the reach of any realistic microscopic 
calculation based on spin-wave theory or other standard condensed matter methods 
methods, such as Schwinger boson mean field theory. The fully systematic effective 
Lagrangian method thus clearly proves to be more efficient than the complicated 
microscopic analysis. Another virtue of the effective Lagrangian technique is that 
it addresses the problem from a unified and model-independent point of view based 
on symmetry - at large wavelengths, the microscopic structure of the system only 
manifests itself in the numerical values of a few coupling constants. Therefore the 
effective Lagrangian method is ideally suited to understand similarities and differences 
in the structure of the low-temperature series for antiferromagnets in three and four 
dimensions, based on symmetry considerations only. 
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A Evaluation of the cateye graph in d=2-\-l 



The singularities contained in the integrals Ji and J2, originating from the cateye 
graph 5c, may be removed by subtracting suitable counterterms. Since our main 
focus is the 0(3) antiferromagnet, here we only discuss the renormalization of the 
function J2 - the renormalization of the quantity Ji, which does not contribute to the 
free energy density for A^=3 according to eq. fl3.14p . will be discussed elsewhere. 

The singularities contained in J2 may be removed by subtracting the following 
counterterms: 

J2 = J2 - - c^gi{M,T). (A.l) 

To establish this result, we use a method, developed in [23|, which at the same time 
also provides us with a representation of the renormalized integrals suitable for nu- 
merical evaluation. We first cut out a sphere S around the origin of radius |iS| < P/2 
and decompose J2 accordingly: 

J2 = jd'^x [d^G{x) d,G{x)y + j^d^x [d^G{x) d,G{x)Y ■ (A.2) 

In the integral over the complement T\S of the sphere, the integrand is not singular 
and the limit c? — 3 can readily be taken. In the integral over the sphere, we insert 
the decomposition (14. ip : 



J2 



= jd'^x (^d^Gd^^oY + ^d^Gd^,Gd,Gd,A + Ad^Gd^AdMd^A 
+ 2d^G d^G d,A d,A + 4d^G d^A d,A d,A + js^A 9^a}' j . 



(A.3) 



In d=2+l the first four terms are convergent. However, the last two terms, involving 
three and four non-thermal propagators, respectively, are divergent. 



In order to extract these two singularities showing up in d=2+l, we follow Ref. |23| . 
where the dimension was d=3+l. We first disregard derivatives and consider the 
expression AGA'^, which contains three non-thermal propagators. Since A{x) is Eu- 
clidean invariant, the integral 

Va;4G'A^ (A.4) 
s 
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only involves the angular average of G{x) 

f{R)= [d'^-^nCix), R=\x\. (A.5) 




The differential equation 

(A.6) 

implies 

= 0, R<p. (A.7) 

Since / is regular at the origin, the differential equation fixes it uniquely up to a 
constant. The function gi ch{Mx4) obeys the same differential equation as G{x) and 
coincides with it at the origin. The angular averages of these two quantities are 
therefore the same, i.e., 

d'^s GA^ = gif d'^x ch{Mxi) . (A.8) 
s J s 

We split the integral over the sphere into two pieces, 

Agi [d'^xch{Mx^)A^ = 4:gi [d'^xch{Mx^)A^ -Ag^ [ d'^x chiMx^) A\ (A.9) 

where the singularity is now contained in the integral over all Euclidean space, in the 
form of the counterterm 

C2 = 4: [d'^xch{Mxi) A^ . (A.IO) 
Jn 

The same line of reasoning goes through for the expression 49^(5(9^ A d^^Ad^A in 
eq. flA.3p . where one ends up with the counterterm 



C4 = 4 / d^xd^chiMx^) d^A d^Ad^A . (A. 11) 

As far as the last term in eq. (lA.3l) . involving four non-thermal propagators, is con- 
cerned, it suffices to subtract the temperature-independent integral of |(9^A(a;)9^A(x)| 
over all Euclidean space, 

C3 = jd^x {d,Ad,Ay , (A.12) 

in order to remove the singularity. Collecting the various pieces, we thus arrive at the 
following representation for the renormalized integral in d=2+l: 



2 



^xT+ / d^xU- / d^-^ 
r Jt\s Jn\s 



J2 = ld'xT+ I d-'xU- I d-'xdf.Adf.A-W , 



T = (^df^Gd^Oy + 4df,Gdf,Gd,Gd,A + 4d^Gdf,Ad,Gd,A 

+2d^Gd^Gd,Ad,A, 
U = 4d^Gd^Ad,Ad,A + d^Ad^Ad,Ad,A, 

W = Agid^ch{Mxi) d^A + d^Ad^A. (A.13) 
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This expression involves ordinary, convergent integrals. Exploiting the fact that G{x) 
and A{x) only depend on r = |x| and on t = X4, the integrals occurring in this 
representation become effectively two-dimensional 

d^x = 27rrdrdt. (A. 14) 

Note that the quantity J2 must be independent of the size of the sphere - this provides 
us with a welcome numerical consistency check of our calculation. 

It is instructive to compare our decomposition of the integrals (1A.13P with the 
decomposition originally used in Ref. 23|], which in d=2+l amounts to 



J2 

U 
V 

w 



d^xU+ /d^xF- 

T\S 

d,Gd,G 



d'x (9^ A ■ W . 



■R.\S 



d^G S^g) ' + Ad^Gd^G d,G d,A + 2Q^^ d,A d,A + AQ^, d^A d,A , 
w + Agid^ch{Mxi) d^A + d^A d^A , (A. 15) 



with 



G 



fx" - ^x'xi + 9xt) + UM'xt go gi + 2(2M"xt + M'^x'xi) gf 

d^G{x) di,G{x) — G^aGyp x^xp , 
-\5^ugo + 5l5t{\go + M^g^). 



(A.16) 



The main difference between the two decompositions ( 1A.13I) and (1A.15I) concerns the 
terms involving two thermal propagators, where we have 



d'x [2d^G d^G d^A d^A + 49^G (9^A d^G d^A 
= I d^x (2Q^^ dyA dyA + AQ^y d^A 9, a) - j d\ d^A d^A 
+ y d=^x (2G^,^G^p x^xp dyA d,A + AG^^G d^AdyA). 



w 



(A. 17) 



Now, in four dimensions the integral over the sphere on the left hand side contains a 
logarithmic singularity - this was the reason why in Ref. 23[ the above decomposition 
was performed: the singularity then occurs again on the right hand side in the integral 
over all Euclidean space. It turns out that, in order to renormalize the integral J2 
in d=3+l, it is thus not sufficient to just subtract the two counterterms C3 and C4 in 
eq. (lA.l[) - rather, one has to subtract two more terms, 

J2 = J2-cs-c,gi + Ud + Q){d - 2) A {G^^f + - 2) XM\g,f , (A.18) 



in order to remove all singularities in J2. 
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In three dimensions, as we have seen, the integral over the sphere on the left hand 
side of eq. flA.lTp is not singular. Likewise the integral over all Euclidean space on the 



right hand side is perfectly well defined, such that there is no need to introduce the 
above decomposition ( 1A.15I) in three dimensions in the first place. Still, in order to 
check that our entire calculation is consistent, we have verified both analytically and 
numerically that the evaluation of the quantity J2 via eqs. (lA.13p and (lA.lSp yields 
the same result. 

We close this section with a comment regarding dimensional regularization and 
the different structure of the singularities in three and four dimensions, respectively. 
The essential point can be seen in the identity ( 1A.17I) which involves two thermal 
propagators. In four dimensions, as we have seen, the cateye graph is of order and 
so is the singularity occurring in the last term of eq. (lA.17p . Now, in four dimensions 
we also have two- loop graphs which are of the same order p^: graphs 8D and BE 
(see Figl2]) which involve vertices from the next-to-leading order effective Lagrangian 
Therefore, the singularities in the last term of eq. (1A.17P - the integral over all 
Euclidean space - can be absorbed into a combination of these next-to-leading order 
coupling constants. 

In three dimensions, on the other hand, these two-loop graphs are of order p^, 
i.e., beyond the order p^ considered in the present work. Loosely speaking, there is 
no communication between the cateye graph (order p^) and these two-loop graphs of 
order p^ and it so seems that - in three dimensions - the 'singularities' in the last 
term of eq. (lA.17l) cannot be absorbed, as there are no next-to-leading order coupling 
constants available. However, in three dimensions the parameter A arising in the 
last term of eq. flA.17p is finite, such that the bookkeeping of 'divergences' in three 
dimensions is perfectly consistent. 



B Renormalization 

In appendix B, we would like to derive the expression fl4.9p for the free energy density 
of an 0{N) antiferromagnet in d=2+l. We first show that all the terms in the free 
energy density (13. 140 that are linear in the kinematical functions gr{M,T) can be 
merged into a single such function, namely go, by replacing the bare mass M with the 
physical mass M.^. 

With the decompositions (14. 2 p and (14. Sp . the terms in the free energy density 
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fl3.14p linear in gr{M,T) read 
^{1} . 



-l{N-l)go{M,T) + l{N-l){N-3) — Xg,{M,T) 



- l(iV - l)(iV - 3)^ — ^72(M, T) + {N- - h) — g,{M, T) 

M"^ 1 
+ ^C2(iV - l)(Ar - 3) _ g,{M, T) - \c,{N - 1){N - 2) — g,{M, T) 

+ i(iV-l)(iV-3)(2iV-5) — AV(M,T). (B.l) 

Now the pressure at low temperatures is of order exp(— M^/T) originating from one 
particle states - states containing two or more Goldstone bosons only show up at 
order exp(— 2M^/T). Therefore it is possible to extract the physical Goldstone boson 
mass from the behavior of the pressure at low temperatures, 



Using the relation 



M^ = -limTlnP. (B.2) 



9r^i = , (B.3) 



this limit amounts to 

Ml = M^ + (N-3)X^ + {2(k, - *0 + |r + '-^} 1^ + 0(M=) , (B.4) 

where the coefficients bi and 62 are given by 

h = ^(iV-3)72-i(iV-2)74, 

62 = (iV-3)(2iV-5). (B.5) 

The quantities 72 and 74 are singular functions of the dimension d and are related to 
the coefficients C2 and C4 - defined in eqs. (lA.10p and (lA.lip - as follows: 

C2 = 72M2'^-6 , C4 = 74M2'^-2 _ (3 

Inspecting the curly bracket in formula ( ]B.4p one thus notices that the infinities con- 
tained in C2 and C4, which stem from the three- loop graph 5c, are absorbed into the 
combination k2 — ki of next-to-leading order coupling constants, originating from the 
one-loop graph 5d. Note that in d=2+l the parameter A is finite, such that second 
term on the right hand side of eq. (lB.4p . coming from the two-loop graph 4b, does not 
contain any singularities. 

One readily verifies that the replacement gQ{M,T) go{MTr,T) in the first term 
of eq.l lB.ip cancels all other terms linear in gr{M,T). The free energy density, linear 
in the kinematical functions, thus takes the simple form 

-i(iV-l)(7o(M^,T). (B.7) 
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We now proceed with the terms in the free energy density that are quadratic in 
the kinematical functions gr{M,T). They are 



z 



{2} 



|(iV - l)(iV - 3) — gr{M, Tf + ^^(iV - l)(iV - 3)(7iV - 17) — A g,{M, Tf 
- \{N-l){N-?,f — \g,{M,T)g2{M,T). (B.8) 
In the first term we make the replacement gi{M,T) gi^M^^jT), which amounts to 

g,iM^, Tf = g,iM, Tf - |2(iV - 3)^ A + 0{M^) } gi{M, T) g^iM, T) . (B.9) 
One notices that this cancels the third term in eq. flB.Sp . We are thus left with 

i^N-l){N-3)—g,{M^,Tf + ^^{N-l){N-3){7N-17)yl-Xg^{M^,Tf. (B.IO) 

Note that in the second term we have also replaced the bare mass by the physical 
mass, both in the prefactor and in the kinematical function: this is legitimate as 
the difference is beyond our accuracy. Finally, in the prefactor of the first term, we 
also express the bare mass by the physical mass using relation (lB.4p . obtaining the 
following terms in the free energy density quadratic in the kinematical functions: 

l^N-l)iN-3)j^g,{M^,Tf + ^^{N-l){N-3)i5N-U)y^\g,{M^,Tf. (B.ll) 



To end up, we take care of the remaining terms in the free energy density that are 
either cubic in the kinematical functions gr{M,T) or are related to integrals over the 
torus: 

- 

+ TsiN 

Again, we replace the bare mass by the physical mass in the above terms, both in the 
kinematical functions and in the prefactors, as the difference is beyond our accuracy. 
No cancellations of terms occur here. 

Collecting the various contributions, we arrive at the expression for the free energy 



-l)(N-3){3N-7)—g,{M,Tf 

-l)(N-3f^g,iM,Tfg,iM,T) 

- 1 - 

-^)iN-3) — .h- \{N -l)[N-2) — .h. (B.12) 
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density of an O(A^) antiferromagnet in d=2+l: 



zo- UN -1)90 + \{N~l){N-^)^{g,f 



- vkiN - l)(iV - 3)(5iV - 11) ^ (^0' + UN - l)(iV - 3)(3iV - 7) 1^ (^^i 



- UN'^){N-?>f^{g,fg2 + ^(iV - l)(iV - 3) Ji 



- i(7V-i)(iV-2)-lj2 + 0(/). (B.13) 

Note that only the physical mass M^r occurs in the above formula: in particular, 
the kinematical functions are = (y'r.(M^,T). Remember that the temperature- 
independent contribution zq is the vacuum energy density given in eq. fl4.6p . 



C Properties of the kinematical functions gr(M, T) 

In this appendix we discuss some properties of the kinematical functions gr{M,T), 
defined by 

gr{M,T) = 2 / -^^/-^ exp(-pM2) J] expi-n'/ApT') . (C.l) 

We follow the appendices A and B of Ref. [l| where the analogous kinematical func- 
tions for finite volume were considered. Here we adapt the method described therein 
to finite temperature. 

We are particularly interested in the expansion of gr{M,T) in the limit M — 0, 
where infrared singularities occur. We introduce the Jacobi theta-function, 

oo 

S{x) = ' (C-2) 

n=— oo 

and express the kinematical functions ( IC.ll) by 

rpd-2r roo / fuf'^f \ r -i 

9r{M,T) = j^l dtf-l~^ exp ^-^j [Sm - l] . (C.3) 

The integration is split into two regions, < t < 1 and 1 < t < oo. In the second 
region we use the identity 

S{x) = ^S{l/x) , (C.4) 
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change the integration variable t 1/t, and arrive at the following representation for 
the kinematical functions gr{M,T): 



9r{M,T) 



T 



•d-2r 



{Airy 1^^ '■-f+s "'■-§1 ' 



(C.5) 



with 



ttr = I dtf 2 exp 
'o 



'47rT2 J . 



S{l/t)-l 



+ I dtt ''+2 2 exp 



' A-kTh) 



and 



"(C.6) 
(C.7) 



The function dr does not contain infrared singularities and the expansion in powers 
of is of the form 



E 

n=0 



M2 1 



where 



dtt^~^^S{l/t)-l^ 



{C.i 



(C.9) 



The infrared singularities are contained in the incomplete F-function hg'- 

\ ^/ N \ ^ 1 



n=0 



47rT2/ n + s' 



(C.IO) 



The pole in the function T{s) at s = 0, —1, —2, ... is compensated by a pole occurring 
in the second piece of hg which is analytic in M. The two singularities can be merged 
and one ends up with a logarithmic contribution. Details can be found in Ref. - 
here we give the explicit expression for h-N (A^ > 0), 



N 



-N 



71=1 



(-1)^+1 f \N ( M2 ^ 1 



y- 



\n 1 

A^J N-n ' 



(C.ll) 



where the quantity 7 ~ 0.577 denotes Euler's constant. Note that, in the second sum 
over n, the value = is to be omitted. Moreover, for A^=0 the sum in the curly 
bracket is to be omitted. 
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We now consider the kinematical function go{M, T) both in three and four dimen- 
sions. According to eq. flC.Sp we have 



gtHM,T) = 
gt\M,T) = 

Using relation ( 1C.4P and the identity 

TT-^/' r(z/2) c(^) = ^ dtr^-' {s{t) - 1} 



(C.12) 



(C.13) 



one readily shows that the various contributions in eq.( ]C.12!) can be merged into a 
single series in M, involving Riemann zeta functions, 



gtHM,T) = 



C(3) 1 M2 r M2 



1 M3 



^ (-1)" / M \ 2" 



at (M.T) 



rp4 



n=2 

45 



(T > M) , 



1 M2 1 M3 (27 - I) 1 M\ M' 



12 T2 67r T3 



n=3 



327r2 T4 327r2 167r2T2 

(C.14) 



Explicitly, the first few terms in these series read 

1 



TC 



1M2 1M2 M2 



6 ^ 96 ^ 



(T>M) 



TC 

45 



15 M2 15 45(7 - f - ln47r) 
1 - — — H H '4 



47r2 T2 27r3 



167r4 



2^4 



+ 



45 , M2 
m 



327r4 T2 



-(f)' 



(C.15) 



The expansion for the function gQ^^{M,T) coincides with the expression derived in 
the appendix of Ref. [23], where a different method was used. The series for the 
kinematical functions gi{M, T) and g2{M, T) can readily be obtained using the relation 



9r+l 



dgr_ 

'dM2 



(C.16) 
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